Reconstruction of Black Hole Metric Perturbations 
from Weyl Curvature II: The Regge- Wheeler gauge 
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Abstract. Perturbation theory of rotating black holes is described in terms of the 
Weyl scalars tp4 and V'o! each satisfying the Teukolsky's complex master wave equation 
with spin s — =f2, and respectively representing outgoing and ingoing radiation. We 
explicitly construct the metric perturbations out of these Weyl scalars in the Regge- 
Wheeler gauge in the nonrotating limit. We propose a generalization of the Regge- 
Wheeler gauge for Kerr background in the Newman-Penrose language, and discuss 
the approach for building up the perturbed spacetime of a rotating black hole. We 
also provide both-way relationships between waveforms defined in the metric and 
curvature approaches in the time domain, also known as the (inverse-) Chandrasekhar 
transformations, generalized to include matter. 



PACS numbers: 04.25.Nx, 04.30.Db, 04.70.Bw 
1. Introduction 

There is a formulation of the perturbation problem derived from the Newman-Penrose 
formalism [201 ^^^^ is vahd for perturbations of rotating black holes. This formulation 
fully exploits the null structure of black holes to decouple the curvature perturbation 
equations into a single master wave equation that, in Boyer-Lindquist coordinates 
(t, r, 9, ip), can be written as: 
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where M is the mass of the black hole, a its angular momentum per unit mass, s the spin 



of the perturbation, E = r^+a^ cos^ 9, and A 



'^—2Mr+a'^. The source term T is built 
up from the energy- momentum tensor ^^j. Gravitational perturbations, corresponding 
to s = ±2, are compactly described in terms of contractions of the Weyl tensor with 
a null tetrad. The components of the tetrad (also given in Ref. [23]) are appropriately 



dt 
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chosen along the repeated principal null directions of the background spacetime [see Eq. 
(121) below]. The resulting (infinitesimal) coordinate and tetrad invariant components of 
the Weyl curvature are given by 

^ ^ ; Pa'V4 = -px'^Cnmnm for s = -2 
^ tpo = -Cimim for s = +2 

where an overbar means complex conjugation and px is given in Eq. (jSJ below. 
Asymptotically, the leading behavior of the field \E' represents either the outgoing 
radiative part of the perturbed Weyl tensor, (s = —2), or the ingoing radiative part, 
{s = +2). 

The components of the Boyer-Lindquist null tetrad for the Kerr background are 
given by 

(Ik") = (3a) 

(^k") = 7yr^-\ ^ {r' + a\ -A, 0, a) , (3b) 

2 [r-^ + a-^ cos"' u) ^ ' 

{m^f^) = —= ^- (msin^, 0, 1, z/sin6') . (3c) 

V2(r + iacosO) 

One can also define directional derivatives as 

6 = m^d^ ■ A = UK^S^ ; i) = \^^d, . (4) 

With the above choice of the tetrad the non-vanishing spin coefficients are (where 
an overbar stands for complex conjugation) 

1 _ cote 

Pk - - 7 : ^, Pk - -Pk—j=, 

[r — ta cost!) 2v2 

2 sin 6* _ sin^ 

ttk = iapK-j=-, tk = -lapKpK-^, 

2 _ A 

Pk = PkPk y 5 aK = T^K- Pk, 
_ jr-M) 

7 =P + PkPk 2 ' ^ 

and the only non- vanishing Weyl scalar in the background is 

^2 = Mpl. (6) 

As we mentioned above the Weyl scalars ipi and ipo allow a direct computation of 
the radiation escaping to infinity and going down the horizon. The time-domain 
formulation is particularly well suited for interfacing with full numerical relativity 
techniques jH [H |SJ |21 E] • There are, besides, other physical phenomena of interest such 
as the self-force on a particle orbiting the hole |TB|, studies of the horizon structure, and 
second order perturbations [Hj, that require the computation of the metric perturbations. 
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Starting from the pioneering work of Chzranowski ^U] there is a series of papers 
|151 n7\ 1211 1^ deahng with the problem of metric reconstruction by the introduction of 
a potential that satisfies the Teukolsky equation, being neither the ip^ or i/jq describing 
to the physical situation under study. The problem of relating the potentials introduced 
to describe the radiation gauges and the physical ip4^ and ipo has been recently studied 
in Ref. ^H]- There, the results are explicitly given for vacuum metric perturbations on 
the Schwarzschild, i.e. nonrotating, black hole background. 

Given the difficulties in obtaining by these method the explicit metric expression 
for perturbations around a Kerr, i.e. rotating black hole, background, we present 
here an alternative approach. In this paper we give explicit formulae for the metric 
reconstruction in the Regge- Wheeler gauge, still for a nonrotating background, but 
allowing for source terms, bearing in mind, for instance, the applications to the radiation 
reaction problem. In the next section we will give explicitly the form ip^ and ipo take in 
terms of metric perturbation, in the Regge- Wheeler gauge, making explicit use to the 
multipole decomposition of the metric. In order to invert these expressions, in Sec. OJwe 
introduce the symmetric and antisymmetric components of the Weyl scalar under the 
discrete parity transformation. With the help of the field equations of General Relativity 
in the Regge- Wheeler gauge (reviewed in [Appendix A ) we succeed in expressing the 
metric perturbations in terms of and ipQ, including matter terms. In the final section 
of the paper we describe how to generalize the first few of these steps to the Kerr 
background case, and speculate about the completion of this program. 



2. Weyl Scalars 

The first step in explicitly constructing the metric perturbations is actually computing 
the inverse relation, that makes use of the definition of the Weyl scalars Q in terms 
of the Weyl tensor. Chrzanowski [lOj, made this computation explicitly relating the 
perturbed Weyl scalars to the metric perturbations 



(A + 71 + 37 - 7)(A + /i + 27 - 2^)h-r 
(A + /i + 37-7)(5-2r + 2a) 



+ ((5 + 3a + /3-r)(A + 27I + 27) 



^(nm) ? ) 



(7) 



and 
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^0 = -<; - a - 3p + W){5 - 2a - 2p + 7i)hn 

+ {D - p - 3e + e){D - p - 2e + 2e)K 
{D - p - 3e + e){5 + 2n ~ 2p) 



+ {6 - a - 3p + W){D - 2p - 2e) 



h 



(Im) 



(8) 



.nn - n^n'^/i^^, /iim = \^m^h^^, etc. 



where /i^ 

Many simplifications are possible in the analysis when the background has spherical 
symmetry. In the Schwarzschild black hole case expressions dZI)-® reduce to 
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(9) 



(10) 



where / = 1 - 2M/r and /' = 2M/r2 . 

The imposition of spherical symmetry also carries the following computational 
advantage: the multipole decomposition of the metric perturbations in terms of spin- 
weighted harmonics _2yem{9) can be performed [221 HHj, and even and odd parity 
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perturbations decouple so they can be considered independently. Below we shall 
decompose all metric perturbations in multipoles with index im (not to be confused 
with the tetrad vectors). 

From Eq. Q and (fTUj) in the Regge- Wheeler gauge {h{"^ = h^J^ = G^™ = = 
iodd)i^em-^^ we get 



1p4 



and 



/ 



E 



im 
X 



im 



16r2 



(£ + 2)! 



2^imi 



im Im 
1 

X 



4/r2 

i 



^ (^-2)! 
(£ + 2)! 



2/2^2 + - /I (^"""^C + / ^""'^M'") } +2Yim. (12) 

These represent our basic equations (real and imaginary parts) that we will use in 
the next section to express the metric perturbations in terms of ip^ and ipQ. 



3. Explicit solution in the Regge- Wheeler gauge 

Two key elements are introduced here in order to complete the inversion of metric 
coefficients from Eqs. ^JandEl The first is the decomposition of the Weyl fields into its 
symmetric and antisymmetric parts with respect to the discrete parity transformation. 
This allows to separate the even and odd parity perturbations from the multipole 
decomposed V'o, and metric perturbations. For the even parity case this allows 
to obtain directly two of the four metric coefficients. In order to obtain the other two, 
we have to resource to the General Relativity field equations, which represent the other 
key element in the inversion process. For the odd parity case, one ends up with first order 
differential relations that can be brought to explicit integrals (previous simplification by 
making us of the odd parity field equations). 



3. 1 . Even parity 

Let us define the symmetric and antisymmetric Weyl scalar fields as [TS] 



± (-)' 



^—l,—m 



(13) 
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where for notational simplicity we dropped the £, m indexes. Thus, given the symmetric 
nature of the even parity metric perturbations, Eqs. (fTT|) and fll2|) take the form 

and 

^o" = 4^ i^t + ^Hi- + Ht) . (15) 

From Eqs. p4|) and (fT3j) we can obtain the im components of the metric 
perturbations as follows 



Hr{r,t) = - — 



4r2 /(£ + 2) 



/ V (^-2)! 



and 



H',"^ir,t) + H',"^ir,t) 




(£ + 2)! 
(£-2)! 



^4 



(16) 



(17) 



We now bring into the play the Hilbert-Einstein's equations in the Regge- Wheeler 
gauge, Eq. ()A.14|) give, 

167rr^ 



i7o^-(r,t)-i7^(r,t) 



V2A(A + 1) 



where -Ff^is a source term given in Table III of Ref. 29^ (See also [Appendix B ) This 



rim 



, and 



The last metric 



rim 



allows to find the metric perturbations for i^Q , Hi 

coefficient has to be found by use of the Hilbert-Einstein equations [i^^™ and H^ — H, 
give a measure of the trace of the even parity sector in the Regge- Wheeler gauge, so it 
does not appear in the Weyl scalars since the Weyl tensor is traceless.] 

Using Eq. ()A.12|) we can solve for drK^"^ and then replace it in Eq. ()A.8|) to find 
j^im terms of the other even parity metric coefficients and source terms 

K{r, tY"" = 

(r - M) fm"'{r,t) , {r-2M)r^H^"'{r,t) 



r^d Hi^ir,t) , M£Hl^ir,t) 



X 

rim 



jjtml 

A 



d Trim I 



A 



+ 



A 

Im I 



(2 r2 - 8 rM + 9 M^) i7^"^(r, t) 

rA (r - 2 M) 
(-r^A + 2 rMX + 3 - 2 rM) H^'^ir, t) 



rX {r-2 M) 
-3r + 7M) fM,{r,t) 



X{r-2M) 
{7 M - 4:r)nrB^"'{r, t) 



+ 



A£m 
^0 



-2M)TTr^^B^"'{r,t) 

VxTix 

V, t)r^ 



VxTTx 



A (r - 2 M) ' 



(19) 
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This form of the metric coefficient K^"^, involves second derivatives of the Weyl 
scalars. One can consider an alternative integral form (on the hypersurface t = constant) 
derived from Eq. ()A.12|1 



K 



dr 



1 2M 
2M -1- — 



Ira 



dt 



, 2M . 
+ —^0 



(20) 



-16n (r - 2M) 



^7r(r 



v/2A(A + 1) VXTl 



2M)^,„ 



It is worth mentioning here that if the source is modeled as a particle (represented 
by a Dirac's Delta) the above metric coefficients are continuous (C^) for headon 
collisions p£j| at the location of the particle. However, for more general orbits they 
do not all appear to be continuous but some of them behave as a Dirac's Delta. For 
instance, one can see that from expression ()18|l . F^"^ is proportional to a Dirac's Delta, 
as given in Table 



11 hence, at least ifg 



or i?!'" have to behave as 5[r — R{t)]. 



3.2. Odd parity 

From Eq. and (fT^ . given the antisymmetric behaviour of the odd parity metric 
coefficients, we get 



C = ^ W [dt - fdr + /'] {^"''^h'^ - f , (21) 



and 



ro = 27^y |:p|y [dt + fdr - f] r'^h'^ + f . (22) 

A linear combination of these previous equations produces 



C + yV^o- = 4^:^1/1:^11 + (fdr - f) if , (23) 



and 



- ^^0- = i^yl^ [fdt^"'''h{- + {fdr - /') ( • (24) 

From Eq. ()A.17j) we can substitute dth^"^ into Eq. leading to the equation 



9rh'r+[(f) h'r = S'r{r,t) = 
-2ir^ /Ii±2)!/ ■ 



27r ir'^Dim 

/v/A(aTT) 



(25) 
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which integrated produces 



h'r = I r S'rir', t) Jl-^dr' + Ci^it) \ , (26) 

(27) 

where Cf™ is an integration constant, that in vacuum and with vanishing ipQ and ipA can 
be taken to vanish |26j . 

Knowing now the form of /if™ we can use Eq ()24|1 to find a differential equation for 



dr-hQ 



Ira 



(r,t) = 



42r2 /(£ + 2) 



(£-2)! 



im 



(28) 



with solution 



C = (l - ^) f^^-' + ^O^-^W} • (29) 

(30) 

Again, Cq™' is an integration constant, that in vacuum and with vanishing ipQ and ip/^ 
can be taken to vanish. 

This essentially completes the work of expressing the metric perturbations in terms 
of the computed ipo and 1^4^ expressed in the time domain. For the even parity case it 
contains second derivatives of the Weyl scalars (unlike the corresponding expressions for 
the radiation gauge [iH]-) For the odd parity case, solutions (|29|l and (j26|) are written 
in an integral form. 

A last observation applies here, since the spin weight of the Weyl scalars are s = ±2 
they do not contain multipole modes i = and i = 1, hence we need to give them 
by directly solving the field equations for the metric coefficients. The Regge- Wheeler 
gauge do not completely allow to determine them, because there are one extra degree of 
freedom for i = 1 and two degrees of freedom for i = 0. Zerilli j^H] has made choices to 
fix this extra freedom that allowed him to solve analytically for the metric coefficients. 
In Ref. I'B^ a different choice was made to make those coefficients continuous in the 
headon collision of extreme mass black holes. Finally in Ref. |T2j the metric coefficients 
for £ = 0,1 have been found in the harmonic gauge, for particles in circular orbits. 



4. Discussion of Kerr perturbations 



A possible generalization of the Regge- Wheeler gauge conditions for spherically 
symmetric backgrounds, but where perturbations are not decomposed into multipoles 
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isIZI 

{smeyhgg-h^^ =0, (31) 

he^ = 0, (32) 

sin 6* dg{sm9hte) + d^ht^ = 0, (33) 

sin 61 dg{sm9hre) + d^phr^, = 0. (34) 

The first equation above leads to the condition G^™ = 0. The second gives then 
(odd)j^hn _ Q rjj_^g other two differential conditions are chosen such that they lead to 

{even)f^£rn ^ q ^ (even) f^im ^ ^^^^ (odd) f^em ^ q {odd)-^ha ^ Q be UnCOUStraiued. 

Now we will consider the generalization of the Regge- Wheeler gauge in the Newman- 
Penrose formalism. In this formalism, the first two Regge- Wheeler conditions, Eqs. 1)311) 
and ()32|) . have a simple generalization 

/imm = m'^m^'/i^,, = 0. (35) 

Note that requiring that the real and imaginary parts vanish contains both conditions. 
Obviously, 

/isffi = m^m^/^M. = 0, (36) 

also holds. Note that conditions ()35p and (jHUj) are invariant under type III (spin-boosts) 
transformations of the background tetrad 

l^A^l, n-^A'V m^e^*®m, rn ^ e-^^^ni. (37) 

This is an important feature, since the Kinnersley choice of the tetrad, with the spin 
coefficient e = is just a simple, but arbitrary way of fixing the spin-boost freedom. 
In contrast the convenient choice of the / and n tetrad vectors along the repeated 
principal null directions of the Kerr background allows to single out wave equations for 
the perturbations of '04 and t/'o- 

To generalize the differential conditions (jHHI) and ()Mj) one can resort to the type 
III transformation properties of the b and 5, as well as spin coefficient operators in the 
Kerr background acting on the metric coefficients /;,(im) and /i(nTn)- The objects 

(5 - 2a)/i(ifi) ^ A^{b - 2a)/i(ifi,), (38) 

and 

{b + ^ A-\b + 2^)/i(,^), (39) 

transform as objects of spin- and boost weight +1 and —1 respectively under type III 
transformations of the background tetrad (jHZI). 

In order to reproduce the differential conditions ()33|) and ()34|) in the Schwarzschild 
limit, one can then require 

3? {{b -2a + aT- 67r)/i(ifi)] = 0, (40) 

and 

3? [{5 + 2^ + cf- rf7r)/i(nm)] = 0, (41) 
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Where 3? is the real part and where following 28j we introduced additional terms 
containing spin coefficients with spin ±1 respectively and boost multiplied by constants 
a, b, c, d to allow for a more general choice of the gauge. These constants can be readily 
chosen to facilitate the metric reconstruction or, in other contexts, to impose further 
symmetries or facilitate the numerical integration of General Relativity field equations, 
etc. It also stresses the ambiguities in generalizing the Regge- Wheeler gauge on the 
Kerr background. 

Note also that conditions ()4()|1 and ()41|1 are invariant under type III transformations 
as well. A crucial role in achieving that was played by the spin- transformation 
properties of the constructed object, allowing invariance of its Real part. 

Independently, we can try to proceed along the lines of the previous sections with 
now a simple mode decomposition of the metric coefficients. For instance 

1 



where we decomposed 



7)™ -I- 

'^(Im) =t 



(Im) 



(42) 



(43) 



We can now replace this decomposition directly into Eqs. ((Zj) and (jH)) for tp^ and 
ipQ or the following more convenient form making use of Ref. |23, Eq. (2.11), where we 
have the following identity 

D - {p+l)t + e + qp-'p {6 - pI3 + qr) 
= [5 - (p + l)P - a + qr + W] (^D - pe + qp^ , (44) 
and the identity derived from it exchanging tetrads 1 —>■ n and m — > m. 
A + (p + 1)7 - 7 - g/i + 71 (6 + pa - qn) 

= \5 + (p + l)a + P - qn -t] (^A + p7 - qp^ . (45) 
Using p = 2 and g = in the above identities allows us to rewrite Eqs. ((Tj) and (jH)) 

as 



^4 = -|(5 + 3a + /3-r)(5 + 2a + 2/5-r)/i„„ 

+ (A + 71 + 37 - 7)(A + 71 + 27 - 2^)hr_ 
- 2(A+7i + 37-7)(5-r + 2a) 

+ 2(6 + 3a+^-T){-jj)] 



(46) 
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and 

^0 = !<( (5_a-3/3 + 7f)(5-2a-2/3 + 7f)/iii 

+ (D - p - 3e + e){D - p - 2e + 2e)h^^ 
- 2{D-p-3e + e){S + n-2p) 

-2(5-a-3/? + 7f)(p)];i(i^)|. (47) 

A choice of the symmetric tetrad [Appendix^ will further simplify the appearance 
of the equations. At this point we impose our gauge condition on /i and h^^^—^ ■ Then, 
paralleling the work done in the nonrotating limit, we could make further progress by 
writing explicitly the Newman- Penrose equations in terms of metric perturbations, for 
finally using these equations to obtain decoupled expressions for some metric coefficients. 
The completion of this program remains an open issue and goes beyond the scope of 
this paper. We leave this for future research. 
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Appendix A. Hilbert-Einstein equations in the Regge- Wheeler gauge 

Because of some misprints in the original Zerilli's paper [221 we reproduce here the 
relevant equations for our discussions (See also Ref. [131 123)- 
The metric perturbations on a Schwarzschild background 

ds^= -(l-'^)dt^ + il-'^rUr^ + r^de'' + smedif^), (A.l) 
can be decomposed into spherical harmonics |29j 





Hl"'{t,r) 








Hi"'{t,r) 


(1-^) 










/if(t,r)| 






r^G^"'{t,r) ^ 






G^'"(t,r) % 


r2 Ik^"" sin^ 9 + G^™ Z^m 



Yem{0, if) 



(A.2) 

for the even parity modes. 
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-h^^{t,r) ^ 



C(t,r)-f^ -hnt,r)-^ hnt,r)^^ 



sin 6d(p 



"■0 K^T ' ) de 

h^^{t,r) '^^^ 



de 



sin 9 



hT{t,r) 



sinOd 

de 



hi'^{t,r] 



sinOd 

de 



hi"'{t,r) sin^% 



/i^™(t,r) sin^% 
-h^^{t,r) sin^% 



Yem{0, (fi) 



for the odd parity modes. Where 51^1^ = 5'^^'^" + h^^,. 



Above we used ZeriUi's notation 



d r d 



92 



d"^ d 1 

Wim= { - cot 6'— - ^^-^ 



and 



^ + sm^cos^- 



We will also introduce ZeriUi's A 



Appendix A.l. Even Parity 

The ZeriUi's (C7a)-(C7g) equations with corrections are 



^ 2MV d'^K^"' 1 
r / dr'^ r 



2M 



d_ 
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-4:V27TiA 
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r-2M 



2 Qij^tm 



1 dHi"^ 1 

H \ 

r dr r{r — 2M) 



r{r - 2M) 

r-M dK^'^ 
r{r - 2M) dr ~ 



K 



(A + l) 



em 



dHi 



em 



r-2M dt 



< Tjem 



+ 



(A + l) 
r(r - 2M) 



(A.3) 

(A.4) 
(A.5) 
(A.6) 

(A.7) 



(A.8) 



(A.9) 



(A.IO) 
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lTx{r - 2M) 



5, 



im : 



(A.12) 
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dtdr 
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2M\ d'^H, 
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2(r - M) dE\ 
r{r- 2M) dt 
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-- 1 



r ) dr 



r + M dH^ 
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j^im 



j^im 



Snr'^Fem 
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(A.13) 
(A.14) 



Appendix A. 2. Odd Parity 



The Einstein equation in the Regge- Wheeler gauge for the odd parity sector are [See 
ZeriUi'spU] equations (C6a)-(C6c)). Note the corrections to the source terms.] 



d'^hi"' 2 dh\ 



£m 



drdt r dt 



+ 



AM 2(A + 1) 



t,£m 



2M 



(i-^)v/(aTT)' 

r dt 
1 



dt"^ drdt 

^ 2M^dhf^ 
r dr 



+ 2A(r-2M) 



£m 



imir 



2M)Q 



Im 



dh^^ 



dt 



+ 



2M 



tm 



v/A(A + l)' 



(A.15) 
(A.16) 
(A. 17) 



where Q^^, Qim and Dirn give the multipole decomposition of the energy- momentum 
tensor (See Table IXT|) . 



Appendix A. 3. Source terms 

Table lAll gives the source terms produced by an orbiting particle in the Schwarzschild 
background after decomposition of the Stress-Energy Tensor into tensor harmonics. 
There, U^{t) = dt/dr, is the zeroth component of the four- velocity of the particle and 
Qp{t) is its angular location. 



Table Al. Energy-momentum-Stress Tensor in terms of Tensor Harmonics 
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AeUr,t) = moU^t) (f)'(r- 


2M)-25[r - R{t)]Yem{^p{t)) 




2 

r-^5[r-R{t)]Y,^{n^{t)) 


A[ll = V2^moU%t)^r- 


-H[r - i^(^)]F,„(^],(^)) 


B^:l = [X + lr^'^moU^t)\ 






Bem = [A + l]-'/^moU\t){r - 2M)-'^6[r - R{t)]^{Q^{t)) 


X 


f l]-^'^m^U\t) 
sine 9$ dt 


9e rft 


r - R{t)] 


Qem = [A - 

X 


f l]-^/'^mof/°(t) 

1 aF,„rfe 

sin e 9$ 


^(r-2MrM 

— sm e— — 

dQ dt 


r - R{t)] 


Den. = -[2A(A + l)]-^/2^mot/°(t)(5[r - R{t)] 


Fem = [2A(A + 
fd<^dQ— 1 


l)]-i/'mof/°(t)5[r - i? 


it)] 





Appendix B. Reconstruction in terms of metric perturbations waveforms 

Here we recall the metric reconstruction in the original Schwarzschild perturbations 
approach based on waveforms for the even and odd parity perturbations (Zerilli's and 
Regge- Wheeler respectively). We first introduce the gauge invariant expressions for 
these waveforms. We then make use of the general relativistic field equations, in the 
Regge- Wheeler gauge, to solve for the metric perturbations, including nonvanishing 
matter terms. 



Appendix B.l. Even Parity 

We consider the following waveform |17j in terms of generic metric perturbations in the 
Regge- Wheeler notation 



^hri _|_ 



+ 



(A + 1 
r - 2M 
Ar + 3M 



2M 



Ar + 3M 



(B.l) 
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This is related to Zerilli'sj^ even parity waveforms V^f^^ by 

riV2r2 {r-2M)A 
(A + l)(Ar + 3M) 



where for an orbiting particle 

41 = rmoV2 Y,^5[r - R{t)l (B.3) 



and it relates to Moncrief'sJ^ waveform 'ip^Mon by 



.im ^ VMou /g 4N 

V^even (A + 1) ' ^"^'^^ 

The tt component of Hilbert-Einstein's equations gives us the Hamiltonian 
constraint. In the Regge- Wheeler gauge {hf^ = = = 0) it is given by Eq. 
()A.8|) . Only two metric coefficients [K^"^ and 7^2™) appear in this equation and none of 
its time derivatives. Considering the Regge- Wheeler gauge, the definition of ip^^^,^ (see 
Eq. ()B.1|) . and the Hamiltonian constraint, Eq. ()A.8|) . we can express these two metric 
coefficients in terms of ipi^,^ (and source terms) only 

_ 6M' + 3MXr + X{X + iy ^(,_^\ .iM 

^ ~ 9 /■ \ , ni\ ^\ reven I '-'r ^even 



(B.5) 



and 



r2(Ar + 3M) 
(A + l)(Ar + 3M)^ 

_ 9M3 + 9AMV + 3A^Mr^ + A^(A + ly 
2 ~ r2(Ar + 3M)2 "^^^^ 

3M^ - XMr + Xr\ , ^...^2,tm 

+ r(Ar + 3M) ^'''^'^^^ + ~ ^^^^'^-n 

_ 8^^' g ,(0) 

(A + l)(Ar + 3M) ^^"^ 

87rr3(AV2 - 2 A + 10 A rM - 9 rM + 27 M^) (q) 
^ (A + l)(r-2M)(Ar + 3M)2 ^ ' 

From Eq. ()A.9jl and the expressions for dtK^^ and dtH^ in terms of dtip^^^, we 
find the if^'" metric coefficient in the Regge- Wheeler gauge 

- rC.(a^e.enj + _ 2M) (Ar + 3M) 

Svrr^ „ (0) 4v^z7rr2 



(A + l)(r-2M)(Ar + 3M) ' ^"^ (A + 1) 
These equations together with 



9tAZ + -T^-^AZ (B.7) 



Tjem _ Trim , IGTTr^ Fi>m ,^ „s 
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give us all metric perturbations on the chosen hypersurface in terms only of V'e^n 
^tV'even (^nd the source). (See also Ref. [H] for general source expressions) 
For the specific case of interest of a pointlike particle we have 



^ 6M' + 3MAr + A(A + l)r^ _^ / _ 2M , 



r2(Ar + 3M) 
InmoYernit) U%t){r-2Mf 
(A + l)(Ar + 3M)r 



r 



5[r - R{t)] 



(B.9) 



Trim 
^2 



+ 
+ 



9M3 + 9AMV + SA^Mr^ + X'^{X + l)r 
r2(Ar + 3M)2 

^(;,^^3^) ^r^crcn + - 2M)a,2^,t, 

8TimoYern{t)U%t)il 



em 



2Af A r \ 2^2 



)[AV2 + 2XMr - 3Mr + 3Af2 



(A + l)(Ar + 3M)2 
8nmoYeUt)U'>{t){r-2M)\,^ 



(A + l)(Ar + 3M) 



-<5'[r - R{t)] . 



■5[r - Rit)] 

(B.IO) 



and 



where 



- rdrWe^^J + _ 2M) (Ar + 3M) 
87rmoF£„(t) f/°(t) rp (Ar + M)^^^ _ ^^^^^ 

6'[r~R{t)]. 



+ 



(A + l)(Ar + 3M) 
^7rmoF,™(t) t/°(t) r(r - 2M) ^, 



(A + l)(Ar + 3M) 
^7rmo(rfF^™/rft)(r - 2M)rU°{t) 
(A + l)(Ar + 3M) 



5[r - Rit)] 



jjhn ^ jjhn ^ igTrr^mo f/°(t) angl(t) 5[r - R{t)], 



angl(t) 



) - sin^ e I ) 

dt J \dt J 



yyim ^ d^ dQ—em 
dt dt ' 



X = 29<^f9e - cot^jr , 



em 



sm ^ / 



(B.ll) 

(B.12) 

(B.13) 
(B.14) 
(B.15) 



Appendix B.2. Odd Parity 

We consider the following waveform in terms of generic metric perturbations in the 
Regge- Wheeler notation 
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A 



I hi-{r,t) 



-dth'r{r,t) 



2r , 



2M 



-Kre. 



(B.16) 



This waveform is related to the Zerilh's.SOj and Moncrief 'spTO] odd parity waveforms 

l _ 2Af ' 



„l,odd „l,odd 

VZer — rMon 



by (See Eq. (prT6|) ) 



I odd 



Sn ir{r - 2M)Q^™ 



to the Cunningham et al ^1] waveform by 



1 t 



Im 
G 



(£ + 2)! 



2A(A + 1)^ 



and to the Weyl scalar \E'2 



(^ + 2)! 



(B.17) 



(B.18) 



(B.19) 



(B.20) 



2 8(£-2)! r3 ■ 

[Here we used the Kinnersley tetrad, in the Schwarzschild background, and decomposed 
^2 into spherical harmonics]. 

One can use the field equations to write the metric perturbation in the Regge- 
Wheeler gauge 



1 



/^r(r,t) = -(l ^ 



2A/f ATrr'^C)'^^^ 



/i{™(r,t) = - 



1 



Av/(aTT) 

ATxir^Qim 



(B.21) 



For a source term represented by a particle the corresponding metric perturbations 
in the Regge- Wheeler gauge are 



C(^,t) = i(i-^)9.Kdd) 



+ 



hf'{r,t) 



47rmor(r - 2M)U°{t)ang{t)5[r - R{t)] 
A(A + 1) 

1 T* 

2 (1 - 2M) *^°dd 

47rmor3f/°(t)(|i?)ang(t)5[r - 



(B.23) 



(r - 2Af)A(A + 1) 



where 



ang(t) 



1 fde 



sin Q \ dt 



■ptm 



\ — ]d^Y (e,<l>)-sine — (9,$), (B.24) 



dt J 



and -R(t), 0, $ define the trajectory of the orbiting particle in spherical coordinates. 
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Appendix C. (Inverse) Chandrasekhar transformations in the time domain 

Chandrasekhar transformations deal with the expressions that relate the waveforms 
in the metric perturbation picture (Regge- Wheeler and Zerilli's) and the Weyl scalars 
that describe the curvature perturbations (Teukolsky's). Chandrasekhar jH] found the 
transformations in the frequency domain. Those expressions can be generalized to the 
time domain to describe local transformations, and take into account matter terms such 
as a particle. Below we give the explicit expressions for the transformations (we drop 
the (im) superscript in the waveforms for the sake of notational simplicity.) 

Appendix C.l. From Waveforms to Weyl scalars 

To obtain the Weyl scalars from the waveforms ()B.1|) and ()B.16|) for even and odd parity 
respectively we simply substitute into Eqs. (fT^ . (fTHj) . (jSI}, and the expressions 
flB.5|) - ()B.8|) . and ()B.21|) - ()B.22|) for the metric coefficients in the Regge- Wheeler gauge. 
The result is 



levrr^ / (0) (o)X 8z(r-2M)v/27r42 (r,t) 



+ 16 



(Ar + 3M)(A + 1) V ' ^™ " '""J A + 1 

TT r (AV2 - 2 A r2 + 10 A rM - 9 rM + 27 M2) 42 (r, t) 



(A + l)(Ar + 3M)^ 



^ Fem{r,t)V27T (r-2M) ^ 
(A+1) 



A(A + 1) V^^^" — A(A + 1) 

48 ZTT (r - 2 Mf Q,^ ^ IGvr (3r-8M) gfj _ ^_ . ^^ 2) 



A v^ATIr A v/AT1 



le^rr^ / .0) , . ,(0)\ , 8z(r-2M)v/27r42(r,t) 



(Ar + 3M)(A + 1) V ' £m ' " «my ■ 

TT r (AV2 - 2 A r2 + 10 A rM - 9 rM + 27 M^) 42 (^, ^) 
(A + 1) (Ar + 3M)^ 
g F,^ (r,t) V27r (r-2M) 
(A+1) 



(C.3) 
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'(£-2)! 



4/2r y (£ + 2)! 
IGnr"^ 



odd 



A(A 



lQinr{r - 2M) 
A(A + 1) 

48z7i {r-2MYQem , 167r (3r - 8 M) _ ^ 



Im I 



+ 



A VA + lr A VA + 1 

where we introduced the Chandrasekhar notation for 



and 



V- 
w- 
w- 

s- 



2 _ 2M^ [A^(A + iy 



r3(Ar + 3M)2 



2(1- 
2 
2 



2M\ /A + 1 



M 



(Ar2 - 3AMr - SM^ 

r2(Ar + 3M) 
(r - 3M) 



^TT (r - 2 M) 



l-{rQ^^){r,t))-ir^^Q{r,t) 



(C.4) 

(C.5) 
(C.6) 
(C.7) 
(C.8) 

(C.9) 



Av/(aTT) 

is the source term for the Regge- Wheeler wave equation. 

These relations are local, and in the time domain. Compare to the expressions in 
the frequency domain of Ref. 0, Eqs. (345) and (353) in Chapter 4. 



Appendix C.2. From Weyl scalars to Waveforms 

To obtain the inverse Chandrasekhar relations we make use of Eqs. fll6p - (Pm) . and (f^ - 
fl26|) for the metric coefficients entering in the definitions of the even and odd parity 
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waveforms, Eqs. ()B.1|) and ()B.16|) . respectively 



and 



So, finally 



+ 
+ 



r^r-2 M) g,Ho (r, t) r^^-^H^ (r, t) 



A (A + 1) 



A (A + 1) 



r (rMA - 3 + r^A + 6 rM) ^Hq (r, t) 

(A + 1) A (Ar + 3M) 
(2 r^A - 5 rMA - 21 M^ + 9 rM) r^gi/i (r, t) 

(A + 1) (Ar + 3M) (-r + 2M)A 
(M^ A r + 2 r^MA^ - 12 r^M - r^MA - 2 r^A + 42 rM^ - r^A^ - 63 M^) Ho (r, t) 



2 (A + 1) (A r + 3 M) (-r + 2 M) A 



-4 



+ 



r^TT V2A + 2 (-5 r^A - 12 rM + 9 rMA + 21 M^) ^25 (r, t) 



(A + 1)'A (Ar + 3M) 



^4^(0) ^) 



+ 4 



A (A + l)(-r + 2M) 
2 



- I)r27r (2 rM - 11 M^ - r^A + 2 rMA) F (r, t) 
(-r + 2M) (A + 1)^A2 
r^-r + 2M)7T-§-^B{r,t) r^Mvr 72^ F (r, t) 



odd 



(A + 1)'/' A 
2(1- f) 



(A + 1) v/A (A + 1)A' 



A 



r (1 — 



ho + So- dthi 



(C.ll) 



odd 



r I 2 M 
A 1 1-^ 



X _ 2M 



9iSi(r',t) Jl-^dr' 
2M V r' 



(C.12) 



(C.IO) 



Appendix D. Symmetric tetrad 



A further algebraic simplification of the expressions can be achieved by choosing the 
background tetrad such that we treat ip^ and ipo on the same footing, thus allowing 
simple linear combinations of the sort ipf ± ipo in the expression in Section El The 
components of the symmetric null tetrad for the Kerr background are given by 



(Is") 

(ns" 
(ms^ 



An, 



/r^ + 

V 2E'"' y2AS 

/r^ + 

V2S'"'y2AS^ 
' ia sin 9, 0, 1, 



, A 



V2{r + 



la cos I 



sin 6 



(D.la) 

(D.lb) 
(D.lc) 
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With the above choice of the tetrad the spin coefficients are 
J^s = , crs =, Xs =, = 0, 



2 sin 6' 



V2' 
_ sin 9 



A 

[M(r2 - cos^ 9) - a^r sin^ 9] 



Ps ^ fJ'S = \l ^Pk 



2V2AE 
75 = 65 — ia cos 9 



3 



2E3 



[(r^ + a^) cos ^ — 2iar sin^ - i / o - \ 
= 2v/2sin» 

where an overbar stands for complex conjugation. 

The Weyl scalars computed with the Symmetric tetrad relate to those computed 
with the Kinnersley tetrad as follows 
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